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We develop the qualitative theory of the solutions of the McKendrick partial 

differential equation of population dynamics. We calculate explicitly the weak 

\j-4 ' solutions of the McKendrick equation and of the Lotka renewal integral equation 

Ph . with time and age dependent birth rate. Mortality modulus is considered age 

Q ' dependent. We show the existence of demography cycles. For a population 

with only one reproductive age class, independently of the stability of the weak 
solutions and after a transient time, the temporal evolution of the number of 
^H' individuals of a population is always modulated by a time periodic function. The 

periodicity of the cycles is equal to the age of the reproductive age class, and a 
population retains the memory from the initial data through the amplitude of 
oscillations. For a population with a continuous distribution of reproductive age 
classes, the amplitude of oscillation is damped. The periodicity of the damped 
cycles is associated with the age of the first reproductive age class. Damping 
^— ^ I increases as the dispersion of the fertility function around the age class with 

^^ . maximal fertility increases. In general, the period of the demography cycles is 

C^ ' associated with the time that a species takes to reach the reproductive maturity. 

^— ^ I Keywords: McKendrick equation, renewal equation, demography cycles, periodic 

solutions, age-structure. 

AMS classification: 92B05, 92D25 



1 Introduction 



K> , The McKendrick equation describes the time evolution of a population structured in 

^ • age. The first time it appeared explicitly in the literature of population dynamics 

5t ! was in 1926 in a paper by McKendrick ^B]. The McKendrick equation is a first order 

hyperbolic partial differential equation, with time and age as independent variables, 
together with a boundary condition that takes into account the births in a population. 
The existence of classical solutions of the McKendrick equation and their asymptotic 
time behaviour is well established, and there exists in the literature of population 
dynamics a large number of surveys. See for example the books of Gushing Pj, 
Webb 1231, lannelfi [Hi, Keyfitz [B], Parkas |H], Kot [El, Gharlesworth (Ij, Metz and 
Diekmann ^^ and Ghu E]. 



Regardless the fact that the existence of classical solutions of the McKendrick 
equation is well established, there is a lack of specific examples and no explicit solu- 
tions of the McKendrick equation are known. This is due to the particular form of 
the boundary condition which is difhcult to handle analytically, ^H] and |H]- 

The McKendrick modelling approach is an attempt to overcome the deficiencies 
shown by the Malthusian or exponential growth law of population dynamics, intro- 
ducing the dependence on age into the mortality and fertility of a population. In its 
simpler form, the McKendrick model does not describe overcrowding effects, depen- 
dence on resources or, in human populations, economic and intraspecific interactions. 
To include these effects, several other models have been introduced and analyzed from 
the mathematical and numerical point of view, [2SI, QHIi ^ and j22| . 

In demography, in order to make predictions about population growth, another 
approach is in general followed. After measuring birth and death rates by age classes 
or cohorts, demographers use the Leslie model ^Hl, a discrete analogue of the McK- 
endrick equation, ^21 and ^^. 

Our purpose here is to find the solutions of the McKendrick equation in the weak 
or distributional sense, to calculate exactly specific examples, and to derive some of 
their properties. To close the gap between theoretical and computational models, we 
compare the stability properties of the McKendrick and the Leslie discrete models, 
unifying both modelling approaches. 

From the point of view of the qualitative theory of partial differential equations, 
we show that the solutions of the McKendrick equation have cycles, (demography or 
Easterlin cycles, |S] and [2]), as observed in the growth of human [2J| and bacterial 
populations |20| . In the Easterlin qualitative approach, ^ and |13| . the period of 
the cycles is estimated to be of the order of the age of two generations, which, for 
human population, is of the order of 50 years. Here, we show that the period of the 
demography cycles is associated with the age of the first reproductive age class, and, 
for human populations, is in the range 10-20 years. The amplitudes of the cycles 
are damped and the damping is associated with the dispersion of the fertility of a 
population around some maximal fertility age. 

In the next section, we review some of the facts about the McKendrick equation 
and we describe the methodology and organization of this paper. 

2 Background and statement of results 

We denote by n{a, t) the density of individuals of a population, where a represents 
age and t is time. Assuming that, within a population, death occurs with an age 
dependent mortality modulus fJ.{a), we have, "^°' ' = —ii{a)n{a,t). As aging is time 
dependent, a = a(i), and is measured within the same scale of time, ■gf = 1, the 
function n(a, t) obeys the first order linear hyperbolic partial differential equation, 

^ + ^ = -,iaMa,t) (2.1) 

where a > and t > 0. To describe births, an age specific fertility distribution 
function by age class b{a) is introduced, and new-borns (individuals with age a = 0) 
at time t are calculated with the boundary condition, 

n(0,t) = / b{a)n{a,t)da. (2.2) 



Supported by data from bacterial and human populations, |2l)| and |14| . b{a) is a 
function with compact support in an interval [a,/3], where a > and /? < oo, and 
fj.{a) is a non-negative function. Equation (|2.1|l together with the boundary condition 
(|2.2(l is the McKendrick equation, ^H]- Knowing the solutions of the McKendrick 
equation, the total population at time t is given by, 

N{t) = / n{a,t)da. (2.3) 



The existence of solutions of the linear equation (|2.1|l with boundary conditions 
(12.21) has been implicitly proved by several authors and goes back to the work of 
McKendrick [TH|, Lotka [n] and Feller 9'. More recently, Gurtin and MacCamy [THI 
proved the existence of solutions of H2.1|l for a class of models where the mortality 
modulus depends on the total population, /i = /i(a,iV(i)). However, due to the par- 
ticular form of the boundary condition (|2.2|) , no explicit solutions of the McKendrick 
equation have been found. 

In the development of the theory of the McKendrick equation, the efforts have 
been on the derivation of existence of solutions, on the conditions for asymptotic ex- 
tinction and for positive equilibrium solutions, as well as, on the existence of stable 
age distributions in finite age intervals, |2S| and [S]- From the point of view of de- 
mography, most of the conclusions derived from the McKendrick model are based on 
numerically constructed solutions from particular initial data. 

In the following, we consider the more general case, 

n(0,i)=/ b{a,t)n{a,t)da (2.4) 



where the fertility modulus b{a,t) is age and time dependent, and for every t > 0, 
b{a,t) has compact support in the interval [a,/3], with < a < (3 < 00. In order to 
simplify the calculations, we sometimes consider that the function 6(a, t) has a natural 
extension as a zero function to the half real line a > 0. This boundary condition is of 
special interest in demography and economy growth models, enabling the analysis of 
the effect of fluctuations of the fertility modulus along time. 

To calculate the general solutions of the McKendrick equation H2.1|l subject to 
the boundary condition H2.4|l . we take the initial data n{a,0) — 0(a), for a > 0. 
In general, 0(0) 7^ J b{a,0)n{a,0)da. Under these conditions, as t increases, the 
boundary condition (|2.2|l or H2.4|l introduces discontinuities in the solutions of H2.1() . 
Here, we assume that 0(a) e L;^^^(R_|_), implying that n{a, t) is also locally integrable. 
In this case, the total population at time i = 0, N{0), is well defined only if 0(a) 
has compact support in R_|_. As we shall see, this does not introduce any technical 
restriction because asymptotic solutions depend on 0(a) in an age interval of finite 
length. In this framework, the Cauchy problem for the McKendrick equation must 
be understood in the sense of distributions — weak Cauchy problem. 

To calculate explicitly the general solutions of the partial differential equation 
(12.1(1 obeying the boundary condition ((2.4|l and prescribed initial data, we use the 
technique of characteristics, ^21 and [H]. Writing equation (|2.1() in the form, J!/. = 
—fj,{a)n{a, t), the solutions of ((2.1|l are also solutions of the system of ordinary differ- 
ential equations, 

dn 

— ^ -fi{a)n 

t _ ^ (2.5) 

dt 



These two equations have sohitions, 

n(a, t) = n{ao, h) cxp i^- f^^ fi{s + qq - to)dsj .^ g, 

a — flo = t — io 

where oq is the age variable at time t = to- The second equation in l|2.6|l is the equation 
of the characteristic curves of the partial differential equation H2.1|l . Introducing the 
second equation in H2.6|l into the first one, we obtain the solution of the McKendrick 
equation for t < a, 

n{a,t) — (p{a — t)exp[— / ii{s + ao)ds 



a — t 



a-i)exp - / fi{s)ds] {t < a) (2.7) 



where cf>(a — t) — n(a — t, 0) is the initial age distribution of the population at the time 
to = 0. For t < a, the solution (|2.7|l does not depend on the boundary condition H2.4|l . 
and if (j){a) G L^(R_|_) or (j){a) € Ljg^{'R+), n{a, t) is integrable or locally integrable in 
a, provided /x(a) G L[q^(R+) and /Lt(a) > 0. 
For t > a, hy (|2.6|l . we have, 



n(a, t) = n(0, t — a) exp ( — / fi{s)ds \ := n(0, t — a)7r(a) (2-8) 

and 7r(a) can be understood as the probability of survival up to age a of the individuals 
of the population. 

In mathematical demography, renewal theories of age-structured populations are 
in general used. Due to the simple structure of the characteristic curves of equation 
(|2.1() . Lotka jTT] has shown that the density of newborns at time t, B{t) :— n{Q,t), 
obeys an integral equation. Introducing H2.7(l and (|2.8|l into the boundary condition 
(|2.4(l . and making the zero extension of b{a, t) to the half real line a > 0, we obtain, 

B{t) = / b(a,t)TT(a)B{t-a)da+ b(a,t)TT(a) \ ~ { da 

b{a,t)B{t - a)da + g{t) (2.9) 



which is the renewal integral equation of demography, first introduced by Lotka |17| 
and developed later by Feller j^j. 

If the solutions of integral equation H2.9|l are known, then, by (|2.7|l and H2.8|l . the 
solutions of the McKendrick equation can be written as, 

r JlisL.A,la-t) (t<a) 
n{a,t)^\ -('^-*)^^ ' ^ ' (2.10) 

\ ■K{a)B{t - a) (i > a) . 

Therefore, once the solutions of the renewal integral equation (|2.9() arc known, the 
solutions of the McKendrick equation (|2.1|) with boundary condition 12.4|l are readily 
derived. Inversely, if we solve the initial data problem for the McKendrick equation, 
the solution of the renewal equation is also easily derived. 



The time independent or equilibrium solutions of the McKendrick equation obey 
the ordinary differential equation ^ = —ii{a)n. Therefore, any equilibrium solution 
has the generic form, 

n{a) = rio exp I — / ^{s)ds 

where no is a constant. Multiplying the equilibrium solution by 6(a) and integrating 
in a, by the boundary condition (|2.2|) . it follows that /^ °° b{a)e~ 'S t^^'^)'^^da = 1. 
Therefore, we define the Lotka growth rate as the number, 

+00 

h{a)e--^<'^^''^'^'^'da. (2.11) 



It is a well known result that the properties of the asymptotic solutions of the McK- 
endrick equation are determined by the Lotka growth number. 

In the following two sections, we calculate explicitly the weak solutions of the 
McKendrick equation (|2.1(l , as well as the solutions of the renewal integral equation 

In sectional we assume that a population has only one fertile age class. This leads 
to the introduction of a new boundary condition n(0,t) = 6i(t)n(a,t), where 61 (i) is 
a positive function and a = a > is the only reproductive age class. This corresponds 
to the choice b{a,t) = bi{t)6{a — a), where 6{-) is the Dirac delta function. In this 
case, the solutions of the McKendrick equation are easily derived and the renewal 
integral equation l|2.9|l reduces to the functional equation, 

Sm^/ ^iW^S^'^(""^) (*<") (2 12) 

^' \ bi{t)Tr{a)B{t - a) (t > a) . ^' ' 

The choice of the boundary condition ri(0, f) = bi{t)n{a,t) has the advantage of 
describing the overall growth patterns of a population as if fecundity were concen- 
trated in one reproductive age class. This enables to discuss effects associated with 
delays in reproduction and, as we shall see, it is important to determine the periods 
of the demography cycles. On the other hand, this boundary condition together with 
the linearity of the McKendrick equation enables to unify the Leslie and McKendrick 
approach, with considerable advantages for the qualitative understanding of growth 
effects in real populations. 

One of the main results of section 13 is that the solution of the McKendrick equa- 
tion is the product of an exponential function in time by a periodic function with a 
period equal to the age of the only reproductive age class. The shape of the peri- 
odic modulation depends on the initial age distribution of the population, and the 
amplitude depends on the mortality modulus. 

In sectional we consider the general case where fertile age classes are distributed 
along an age interval, and we derive explicit formulas for the weak solutions of the 
McKendrick and the renewal equations. The asymptotic behaviour of the solutions 
is derived as a function of the Lotka growth rate. Then, we explicitly calculate the 
solutions of the renewal and of the McKendrick equation for the case where b{a, t) is 
constant in the age interval [a, /?]. 

We show that the stability and instability of the solutions of the McKendrick 
equation in the weak sense are determined by the Lotka growth rate. In the present 
theory of the McKendrick equation, the stability and instability of the solutions are 
given implicitly by the real roots of the characteristic equation associated with the 



Laplace transform of (|2.1|l . Pj, |3] and ^T]. An algorithmic procedure to determine the 
roots of the characteristic equation has been recently derived by Farkas [7]. However, 
it is not always possible to locate these roots, [7j- With the approach developed 
here, the stability and instability of solutions of the McKendrick equation are directly 
calculated from the fertility function b{a) and the mortality modulus iJ.{a), parameters 
directly measured in demography, |13j and |14j . 

The technique developed in section|31enables to relate and calibrate the parameters 
of both the McKendrick continuum model and the Leslie discrete model, ^Hjj used 
in demography studies, [14j. It is shown that when the number of age classes of the 
Leslie model goes to infinity, the inherent net reproductive number of a population 
converges to the Lotka growth rate. 

In general, we conclude that if the time intervals in the time series of a population 
equals the age of maximal fertility, this time series has an exponential or Malthusian 
growth. For smaller time intervals, the time series is modulated by a periodic function 
in time. This proves that, in a time scale of the order of 10-20 years, there exists 
cycling behaviour in the pattern of growth of a population. Approximating the age 
distribution of a population by the sum of Dirac delta functions concentrated at 
consecutive age classes, a modelling possibility from the computational point of view, 
the modulation of the pattern of growth of a populations becomes almost periodic in 
time (section EJ . 

At the end of each section, we discuss the qualitative aspects of the solutions we 
have analysed. The main conclusions of the paper for demography and population 
dynamics are discussed in the concluding section [T] 

3 Populations with one fertile age class 

We suppose that births occur at some fixed age a = a > 0, and we assume that the 
fertility function is b{a,t) = bi(t)S{a — a), where S{-) is the Dirac delta function, and 
bi{t) is a differentiable, positive and bounded function of time. For this particular 
case, the boundary condition (|2.4|l becomes, 

n{0,t)^bi{t)n{a,t) (3.1) 

with a > 0. Note that, if a = 0, the Cauchy problem for the McKendrick equation is 
only determined for t < a. For t > a, hy l|3.1|) . ri(0, t) becomes undetermined, as well 
as n{a,t) with t > a. 

To extend the solution H2.7I) of the McKendrick equation as a function of the initial 
data to all the domain of the independent variables a and t, taking into account the 
boundary condition (|3.1|l . we first consider the case t — a. By (|2.8|l . and for t = a, 



n{a,t) ~ n(0,0)exp I — / ii{s)ds 

As newborns at time t — are calculated from the boundary condition H3.1() , to make 
this solution dependent on the initial data, we must have. 



n(a,t) = 6i(0)(/)(a) exp ( — / ^{s)ds ] , if t = a (3.2) 

where n(a,0) = 0(a), see Figure ITTI 



We now introduce the sets Tm — {{a,t) : t < a+ma,t > a+{m—l)a,a >0,t> 0}, 
where m is a positive integer, and Tq — {{a,t) : t < a,a > 0,t > 0}, Figure |0] We 
denote by T„i the interior of the sets Tm- So, given the point with coordinates (a, t), 
we have, (a, t) G T„i, with ?Ti > 1, if and only if, i > a and [1 + (i — a)/a\ — to, where 
we have used the notation [x] for the integer part of x. A point (a, t) £ Tq, if and 
only if, t < a. Note that, the domain of the solution (|2.7|) is the interior of the region 
labelled Tq in Figure O 

For t > a and t < a + a, we take the point (a*, t*) on the line t — t*. Therefore, 
t* > a* and t* < a* + a. This point is in the region labelled Ti in Figure 13.11 By 
(|2.6|l , the characteristic line that passes by (a* , t* ) crosses the line a = at some time 
t = t^, and, 

n{a*,t*) = n(0, t^) cxp ( - / fi{s - tl)ds ] 



where tl = t* — a*. Imposing the boundary condition H3.1|l on this solution, we obtain, 
n{a*,t*) = bi{tl)n{a,tl) exp [- J*, /i(s - tl)dsj. By hypothesis, as tl = {t* ~ a*) < 
a, we are in the conditions of solution H2.7|l for t = tl and a = a, and we have, 

n{a*,t*) — bi{t* ~~ a*)n{a,t* — a*)exp I — / fj,{s — tl)ds \ 



= hi{t* ^ a*)(j){a + a* - t*) 

Cnt'-a* pt* \ 

— / /i(s + ao)(is — / ii{s + a* ~t*)ds 

Jo Jf-a* J 

= &i(i*~a*)0(a + a*-i*) 



X exp — / [i{s + a + a* — i*)ds 




(3.3) 



X exp 



Therefore, we have shown that, for t > a and i < a + a, 
n(a, t) — 6i (t — a)(j){a + a — t) 

X exp — / ii{s + a + a — t)ds — I ii{s)ds 

Note that, for t < a, the solution n{a,t) given by (|2.7|l . (|3.2|l and H3.3|l is in general 
discontinuous when the line t — a is crossed transversally. 

We now proceed by induction. Suppose that, up to some integer k > I, the 
solutions of the McKendrick equation l|2.1fl with the boundary condition H3.1|) can be 
written in the form, 



n(a, t) = (j)[ma + a — t) \\bi{t — a — [i — \)a) 




1=1 

f — a— (rn— l)a \ 

li{s + TOa + a - t)ds (3-4) 



[m — \) I fj,{s)ds — I fi{s)ds 




Figure 3.1: Characteristic curves a — ao = t — to for the McKendrick equation p. 1(1 . 
The interior points of the sets Tm — {{a, t) : t < a+ma, t > a+{m~l)a, a > 0,t > 0}, 
with m > 1, and Tq ~ {{a,t) : t < a, a > , t > 0} are denoted by Tm, with to > 0. 
The vertical hne a = a represents the age of the unique reproductive age class of the 
population. The boundary condition at time t is calculated according to the value 
of n{a,t). At i = 0, n(a,0) = (/>(a). Given an arbitrary point (a*,t*) in the domain 
of the partial differential equation 1(2. l|l . n{a*,t*) is obtained following the solution 
n{a,t) along the dotted line until i = 0. 

where m ~ [{t — a) /a + 1], (a, t) 6 T^, m < k and m > 1. As we have shown, ((3.4(1 is 
true for TO = 1. Suppose now that {a*,t*) e Tk+i- Then, by ((2.t)() . the characteristic 
curve that passes by (a*, t*) crosses the line a = at some time t^ = t* — a* , and. 



n{a*,t*) = n(0, tj) exp - / ^(s - tl)ds . 
Imposing the boundary condition 1(3.1(1 on this solution, we obtain, 

n{a*,t*) — bi{tl)n{a,tl)ex.p i — / fi{s~tl)ds\ . 
As [{t* -a* - a)/a + 1] = [{t* - a*)/a] = k, we have, {a,tl) e Tfc, and by l|T^ . 
n{a*,t*) = bi(tl)n{a,tl)exp \ — / fi{s — tl)ds \ 

( r' \ 

= bi{t* — a*)n{a, t* — a*) exp — / ij{s)ds 

= i f[ bi{t -a-ia)\ (f>{{k + l)a + a* - t*) exp (-S) 



\i=0 



where, 



t* —a* —ka 

/i(s + {k + l)a + a* - t*)ds 



+ k fi{s)ds + / ^{s)ds . 
Jo Jo 



Hence, (|3.4|l remains true for m — k + 1, and we have: 

Proposition 3.1. Let n(a,0) — (l){a) G C^(R+) &e an initial condition for the McK- 
endrick partial differential equation 12. 1\} . with a > 0, t > 0, and boundary condition 
iy.l\) . Assume that nia) e C*'(R+) is a non negative function and bi(t) G C^iTL^) 
is positive. Then, in the interior of the sets Tm, with m > 0, the solution of the 
McKendrick equation Y2.1\) is differentiahle in a and t and is given by: 

a) If {a,t) e To, 

n{a, t) — (f){a — t) exp I — / ^{s + ao)ds 

where oq — a ~ t. 

b) If, (a, t) G T„i and in > 1, 



i{a,t) = Wbi{t - a ~ [i - l)c 



i=l 



l>t — a — {m,— l)a ^ 

X(j){ma + a — t) exp I — / ii{s + ao)ds 







X 



exp I — (m — 1) / id{s)ds — / fi{s)ds) 



where oq = ma + a — t, m — [(t — a)/a + I], [x] .stands for the integer part of x, and 
a > 0. 

Proof As 0(a) e CHR+), /^(a) G C°(R+) and bi{t) G C^{R+), it is straightfor- 
wardly checked by differentiation that, in the interior of the sets Tm, with to > 0, the 
solution obtained by the method of characteristics obeys the McKendrick equation 
(EH). D 

In the construction preceding Proposition l3.1l we have shown that the solutions of 
the McKendrick equation hold formally at the boundary of the sets Tm, to > 0, where 
n{a, t) is discontinuous because, in general, (j){Q) ^ 6i(O)0(q:). To extend the solution 
of the equation 1)2. l|l as stated in Proposition l3.1l to all the domain of the independent 
variables, we introduce the concept of weak solution in the sense of distributions, jl2| . 
We consider the space of test functions that is, the space of functions of compact 
support in R^ with derivatives of all the orders. Let 2?(R^) be the space of test 
functions. If f{x) is a locally integrable function, then f[^] = J ijj{x)f{x)dx is a 
continuous functional in the sense of the distribution in P(R'^), "1^. Taking equation 
(12.1(1 and making the inner product with a function tl^{a, t) G P(R^), as ?/; has compact 
support in R^, and after integrating by parts, we naturally arrive at the following 
definition: 



10 



Definition 3.2. A locally integrahle function n{a,t) (g ^^^^^(R^)) is a weak solution 
in the sense of the distributions of the McKendrick partial differential equation 1^2. 1\) . 

f d d \ 

-—ipia, t) + -^—ipia, t) — fi(a)ip(a, t) n{a, t)dadt — 



for anyip{a,t) £ ^(R^.). 

In the conditions of Definition 13. 21 we have: 

Theorem 3.3. Let n{a, 0) = 0(a) (g L\^^{TL+)) be a locally integrable initial condition 
for the McKendrick partial differential equation 1^2. Ip . with a > 0, t > 0, and boundary 
condition fS.lp . Assume that ^{a) G L;^^^(R_|_) n C''(R+) is a non negative function 
andbi{t) G C"'^(R+) is positive. Then, the weak solutions of the McKendrick equation 
HjP are: 

a) If {a,t) G To, 

n{a, t) = 0(a — t) exp — j ii{s + ao)ds 

where ao = a — t. 

b) If (a, t) G Tm and m > 1, 







n{a,t) = \\bi{t — a — {i — l)a) 



i=l 



f>t — a—{ni—l)a ^ 

X(f){ma + a — t) exp I — / ^{s + ai^)ds 







X exp — (tti — 1) / ^{s)ds — / ^(s)ds 

where oq = ma + a — t, a > 0, m — [(t — a) /a + 1], and [x] stands for the integer 
part of X. 

Proof. We have shown previously that the solutions of the McKendrick equation as 
in Proposition 13.11 hold at the boundary of the sets T,„, with m > 1. However, 
at these boundary points the solutions in Proposition 13.11 are not differentiable and 
must be understood as weak solution in the sense of distributions. So, for n{a, t) as 
in Proposition 13.11 to be a solution of the McKendrick equation in all the domain of 
the independent variables, it must satisfies the conditions in Definition 13.21 for any 
ip{a,t) G P(R^). As by hypothesis 0(a) and ii{a) are locally integrable, we are in 
the conditions of Definition 13.21 and we can verify if n{a,t) given by ProDOsition l3.1l 
can be extended as a weak solution in all the domain of a and t. For that, we write 
the solution in a) in Proposition 13 . II in the form, n{a,t) = $(a — t)e^^(''^*'*\ where 
T{a — t,t) — T(ao, t) is a generic function reflecting the functional dependency of the 
exponential term, and we calculate the integral, 

where $(?/) = 0, for y < 0, ^{y) = (p{y), for y > 0, and <i>(y) G Lj^^CR). Introducing 
the new coordinates y — a — t and x — t, we have. 



RXR4 



■^ - /i(a; + y)?/; ) $(y)e '^^^'^^dxdy . 
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As, T(a — t,t) = T(y,x) — Jq fJ.{s + y)ds, and ^ = ii{x + y) — fJ'{a), almost 
everywhere, the above integral evaluates to. 



d 

RxR 



— Ue-'^'^y^^'A ^{y)dxdy = / ^e"'^^^'"'^ ^{y)dy = 



Therefore, the solution a) is a weak solution of the McKendrick equation. For the 
case b), we introduce the new coordinates y — t — a and x = t. Writing the solution 
in b) in the form, n{a,t) = ^{t — a)e~^*^*^'^'*\ and as m = m{y), it follows that 
^ = iJ,{x — y) = n{a)- By an analogous calculation, we have / = 0, and the solution 
b) is also a weak solution of the McKendrick equation. D 

For the particular case a = 0, Theorem l3.3l gives the explicit solution of the renewal 
functional equation 12.1211 . 

Corollary 3.4. Let n{a,0) — 0(a) (g Lj^^(R^)) be a locally integrable initial condi- 
tion for the McKendrick partial differential equation Y2.1]) . with a > 0, t > 0. Assume 
that ii{a) G L\^^{R+) n C''(R+) is a non negative function and bi{t) G C^(R+) is 
positive. Then the solution of the renewal functional equation ^2.1^) is, 

m 

B{t) = <^(ma -t) J|6i(t- (i- l)a) 

i=l 

X exp — / ^i{s + ma — t)ds 

X exp I —{rn — 1) / fi{s)ds 







where a > 0, m = [t/a + 1], and [x] stands for the integer part of x. 

In the particular case where bi(t) is a constant, bi{t) = bi, we have, in Theorem 
ESand CoroUaryEl 6i(t - a) . ..bi{t - a - {m ~ l)a) = 65". If n{a) ee ^ > and 
bi{t) = bi > are constant functions, the solution of the McKendrick equation is, 

where m, — [{t — a) /a + 1]. For the same case, the solution of the renewal equation 
(ITHll is, 

B{t) = bl(l){sa-t)e-''' 

where s — [t/a + 1]. 

To analyze the asymptotic behaviour in time of the solution of the McKendrick 
equation with time independent fertility function, we define the function s{a,t) by, 

e{a, t) = {t- a) /a + 1 - m = (t - a) /a + 1 - [(t - a) /a + 1] (3.5) 

where [(t — a) /a + 1] ~ m, and ?7i > 1 is an integer. For fixed a and with t — a >Q, 
the function e{a,t) takes values in [0, 1), is piecewise linear and is time periodic with 
period t = a. Then, we have: 
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Theorem 3.5. Suppose that (f>{a) is positive and bounded in the interval (0, a] and 
the fertility function bi is a positive constant. Then, in the conditions of Theorem 
\3.Sl we have: 

a) If, Inbi = /„ ^{s)ds, then, for fixed a and t > a, n{a,t) is time periodic with 
period t — a. 

b) If, In hi > L ii{s) ds, then, for fixed a and t > a, n{a, t) ^ oo, as t —>■ oo. 

c) If, Inbi < L fi{s) ds, then, for fixed a and t > a, n{a, i) —> 0, as t —> oo. Moreover, 
for fixed a, the asymptotic behaviour in time of n{a,t) depends on the initial condition 
4>{a) with a in the interval (0,a]. 



Proof. By (I3.5|l . with oq 
6™ exp 
X exp 
= exp 

X exp 
:= exp 

X(a,i) 



ma + a — t, 

t — a— (m— 1) 



/i(s + ao)(is — (to — 1) / ii{s)ds 
Jo I 



^{s)ds 



TO,(ln6i — / ij{s) ds) j exp I — / fi{s)ds 



t{a,t) 



li{s + a — ae{a,t))ds + / ^{s)ds 



TO,(ln6i — / ii{s) ds) \ exp I — / 



where, 



exp 



exp 



ae{a.t) 



-/ ^i{s)dsjxia,t) 

/x(s + a — ae{a, t))ds 1 
li{s)ds 



(3.6) 



a — Q£(a,t) 



Due to the periodicity of e{a,t) in t, for fixed a and t > a, xio-it) is also periodic in 
t, and we can write the solution b) of Theorem 13. 31 in the form, 



n(a,t) — (f>{a — ae{a,t))x{a,t) 

X exp I TO,(ln6i — / ij{s) ds) ) exp ( — / ij{s) ds 



(3.7) 



Therefore, if In 6i — /„ /i(s) ds = and t > a, the solution of the McKendrick partial 
differential equation becomes oscillatory in time. If, ln6i > L ^{s)ds, then, in the 
limit t —>■ oo, m —* oo, and the population density goes to infinity. If, In&i < 
/g ^j,{s)ds, then, in the limit t -^ oo, the population density goes to zero. By H3.7() . 
the asymptotic distribution of a population depends on 4>{a) with a G (0, a]. D 

If we fix an arbitrary large value of a, and if In 6i > J„ fj.(s) ds, then n{a, t) — > oo, 
as i — > oo, implying that asymptotically there is not a limit in life expectancy. 

Due to the periodicity in time of e{a,t), the functions xio-^i) s-nd 4'{^ ~ Oie{a,t)) 
are time periodic with period t = a, the age of the only reproductive age class. 
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Defining the Lotka growth rate of a population as r = n{a,t + a)/n(a,t), by 
Theorem 13.31 and as [{t + a ~ a) /a + 1] — m + 1, the Lotka growth rate associated 
with the McKendrick equation is constant and is given by, 

r = r — = exp In 6i — / fi(s) ds ] — bi exp — / fj.(s) ds ] (3.8) 

n(a,t) \ Jo ) \ h J 

which coincides with (|2.11|) for the choice h{a) = hi6{a — a). 

The Lotka growth rate (|3.8|l has a simple interpretation. By Theorem 13.31 and 
(|3.8() . for i > a, we have •n(a, t + sa) — n{a, t)r^ , where s is an integer. With t = a 
and a+sa — r, we obtain, n{a, t) = n{a, a)r^'^~°'''" . Defining the Malthusian density 
growth function as nM{a., t) = n{a, a)r(*~"''", and integrating riAf (a, t) in a, the total 
population varies in time according to, 

NM{t) = A^(a)r(*-")/" (3.9) 

which is the Malthusian growth function associated to the McKendrick equation. 
Therefore, within an observation time step equal to the age of the only reproductive 
age class, the solutions of the McKendrick equation grow exponentially in time, with 
the Lotka growth rate H3.8(l . 

The qualitative difference between the asymptotic behaviour of the solutions of 
the McKendrick equation and the simple Malthusian exponential growth without age 
structure is associated with the existence of a periodic modulation in the growth of 
populations. This periodic modulation with period a corresponds to the demography 
cycles observed in real populations, [201 ^^'^ |21j . 

The stability condition for the persistence of non-zero solutions as stated in The- 
orem |^21 is determined by the growth rate r: stability or periodicity of solutions if 
r = I; exponential growth if r > 1, and population extinction if r < 1, provided </>(a) 
is not identically zero in (0, a]. 

Theorem 3.6. If4i{a) is bounded in the interval {0,a\ and the fertility function hi is 
a positive constant, then, in the conditions of Theorem \3.!A for fixed a, and t > a, 

n(a,t) . , , , , , / /"" / N , 

— -z— = X(a, t)0(a - ae(a, t)) exp - / ^(s) ds 



is a time periodic function with period a, and m = [{t — a) /a + 1]. 

Proof The theorem follows by (ESJ, (ESI, (ETj) and (gSJ). D 

By the above Theorem, the population age density and in the total population 
has a periodic modulation in time. These oscillations occur around the Malthusian 
growth curve H3.9|l . 

The amplitude of oscillations of the function n{a, t)/r™ can be easily determined. 
From Theorem 13 . 61 the amplitude of oscillations is, 



A = maXa^^Q^a]x{a.,t)cl){a-ae{a,t))exp{~ J^^l{s)ds) 

- mma(z(o,a] x{a, O0(a - "^(ai 0) exp (- /^ fJ.{s) ds) . 

Assuming that the initial distribution of the population is constant in the interval 
(0,a], and as e{a,t) takes values in [0,1), n(a,t)/r"^ varies in the interval, 

(kg- lo ^(«) ''«e~ ^'' ^^^^ ''^ (t)e~ ^^ '^'^•' ''^ 
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Figure 3.2: a) Time evolution of the solution of the McKendrick equation (|2.1() for 
age classes a = 8 and a — 35, in a population with one reproductive age class a — 25. 
The mortality modulus is /z(a) = 0.05 + 0.001a and 6i = 4.77. b) Distribution of the 
density of individuals as a function of age for t = 10 and i = 100. In both cases, 
we have the stability condition 5i = e^", implying that the Lotka growth rate is 
r = 1. All the solutions have been calculated from Theorem 13. 31 with the initial data 
condition <j>{a) = 2, for a < 100, and 4>{a) = 0, for a > 100. 



Then, for a uniform initial age density of individuals, the periodic modulation has the 
age dependent amplitude, 

A{a) = (t)exp( ^(s)ds j n -exp (- / n{s)ds\\. (3.11) 

In the particular case where /i(a) = /i, the amplitude simplifies further and is given 
by A{a) — 0e^'^"~"'(l — e~'^"). Therefore, the newborns age-class has the maximum 
amplitude of oscillations, ^(0) = 0e^"(l — e~^"), and the amplitude of oscillations 
increases (resp. decreases) when /i increases (resp. decreases). 

Another consequence of Theorem 13.61 is that, in the limit t —> oo, the population 
retains the memory from the initial data through the amplitude of oscillations of the 
growth cycles, (see the discussion in lannelli ^J, pp. 37). 

In Figure 1^21 we depict the time and age evolution of the density n{a,t) of a 
population from a uniform initial age distribution with a maximal age class, and 
Lotka growth rate r = 1. We have chosen the age-dependent mortality modulus 
/i(a) = /io + Mi*^! ^-iid as initial condition the density function with compact support, 
0(a) — 2, for a < 100, and (f>{a) = 0, for a > 100. By (|2.3(l . this corresponds to 
an initial population N{0) = 200. The birth constant 6i has been chosen in such a 
way that r = 1. After the transient time a, and for fixed age, the population density 
becomes periodic in time with period t = a. In Figure 13.2b . the decrease of the 
amplitude of oscillations with increasing ^{a) is in agreement with (|3.11() . 

In Figure 13.31 we show the total population iV as a function of time calculated 
from (|2.3|) and Theorem 13.31 In Figure rOk ) . the Lotka growth rate is r = 1. In 
Figure rOb ). 6i > exp(/p ^{s)ds), and the Lotka growth rate is r = 1.09. In both 
cases, we compare the solution N{t) calculated from Theorem 13 . 31 and the Malthusian 
growth curve H3.9|l . 

If the fertility function bi as considered in Theorem l3.5l depends on time, in general, 
the stability can not be decided in finite time. This follows by a similar analysis to 
the one in the proof of Theorem 13. 51 
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Figure 3.3: Total population number TV as a function of time calculated from Theorem 
13.31 and (|2.3|) . The age of the only reproductive age class is a = 25 and the initial 
conditions are as in Figure 2. In a), ln6i — J^ ^{s)ds = 4.77, where ^{a) = 0.05 + 
0.001a, &i = 4.77 and the Lotka growth rate or stability condition is r = 1. In b), 
bi — 5.2, the Lotka growth rate is r = 1.09 and the number of individuals of the 
population goes to infinity. We also depict the Malthusian growth function l|3.9(l . 
N]\j{t), measured in the time scale of the unique reproductive age class. In all the 
cases and after a transient time, the growth curve of the population is modulated by 
a periodic function with period t = a. 

The main conclusion about the solutions of the McKendrick equations for the 
boundary condition H3.1|l with bi{t) = hi is that the pattern of growth of a popula- 
tion is modulated by a periodic function with a period equal to the age of the only 
reproductive age class. The shape of the periodic modulation depends on the initial 
age distribution of the population, and the amplitude depends on the mortality mod- 
ulus. In time steps of the order of the age of the first reproductive class, the growth is 
Malthusian. Pure Malthusian growth is obtained if the mortality modulus approaches 
zero and the population has an uniform initial age distribution. 

4 Populations with several fertile age-classes 

We consider now that the fertile ages of the individuals of a population are distributed 
in some age interval [a,/3], with < a < (3 < cxd. The boundary condition is now. 



n(0,t) 



b (a, t) n (a, t) da . 



(4.1) 



The constants a and f3 represent the ages of the first and the last reproductive age 
classes, respectively. 

As in (|2.7|) . if, t < a, the solution of the McKendrick equation (|2.1|) is given by, 



n{a,t) = (b{a~t)e-■^o^^'^-+^~t)ds 



(4.2) 



where (/>(a) = n(a, 0) is the initial age distribution of the population. 

The general solution of the McKendrick equation in all the domain of the inde- 
pendent variables is obtained in the following way: 

Theorem 4.1. Let n{a, 0) = 0(a) (g L;^^^(R+)) be a locally integrable initial condition 
for the McKendrick partial differential equation ^2.1}) . with a > 0, i > 0, and boundary 
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condition J^. j| ). Assume that fJ,{a) G Lj^^(R^) n C*'(R+) is a non negative function, 
and that, for every t > Q, b{a, t) is locally integrable and differentiable. Then, in the 
strip S ~ \{fl-,t) : a > 0, < i < a}, the weak solution of the McKendrick equation 
\2.1\) with initial data (j){a) and boundary condition \4-l^ is: 

n.(a,t)J^(«-r''""^''"' '^ ^<" (4.3) 

I 7r(a) j ^ b {c, t — a) (f> {c — t + a) dc , if,t>a, and t < a 

where, 

7r(a) = e^.''iX^)'^" 
b{a,t) = 6(a,t)7r(a) 

and the function ns {a, t) as in i4.S{ ) is locally integrable (ns{a,t) G L\^^[S)). Writing 
ns{a,t) as ns{a,t;(j)), and defining the functions, 

0i+i(a) = nsia,a;(j)i) 

where i >0 and 4io{a) — (j){a), we can construct recursively the general solution of the 
McKendrick equation as n(a, t) = ns{a, t — qa; 4>p), where q = [t/a], and [t/a] stands 
for the integer part of (t/a) . 

Proof. To construct the solution of the McKendrick equation in the strip S, the case 
(a = 0, t = 0) follows from the boundary condition (|4.1() . The case t < a and t > has 
been proved in 1)4.2(1 . We now construct the solution for (a,t) S Ti, Figure mi which 
implies that a > 0, t > a and t < a + a. Let {a* , t*) be the point such that, t* > a* 
and t* < a* +a. By 12.6|l . the characteristic line that passes by {a*,t*) crosses the line 

a = at some time i = t^, and n(a*,i*) = n(0,i*)e ■'*?''* ^ *, where i* = t* - a*. 
Imposing the boundary condition 1(4. l|l on this solution, we obtain, 

n(a*,r)=e~A ''^""'i^'*" / b {a,tl) n{a,tl) da . 



By hypothesis, as tl — t* — a* < a, we are in the conditions of solution ((4.2(1 . and we 
have, 

n(a*,r) = e~A ''^'"*^^''' / b {a,tl) n{a,tl) da 



" b (a, tl) cp{a- tX) e-^°^ M^+a-*t)<is-/,': Ks-tl)ds ^^ _ 
Therefore, we have shown that, for t> a and i < a + a, that is, for (a, t) G Ti, 
n{a,t) = / b{c,t~a)(j){c-t + a)e''^o ^ K^+c-t+a)ds-Jl^,.(s-t+a)ds^^ 

J a 

13 

b{c,t-a)ct,{c~t + a)e- ■^-*+» f'^'^'^'-Io ^(^)'i«dc 
-lo^^^i-^ds / &(c,i-a)^(c-i + a)e^''"'^°'^(")''"dc (4.4) 
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where, 

For a = and t = 



b{c, t-a)=b{c,t-a)e- ■'''' "^'^'^^ 



Ti{0,a)— / b{c,a)n{c,a) dc. (4-5) 

J a 

As, for c > a, n{c,a) G Tqi we introduce l|4.2|l into (|4.5(l . and we obtain H4.;-{|l . 
Therefore, we have constructed the solution of the McKcndrick equation in the strip 
S. 

As 4>(a) is locally integrable, the integrals in 14.3|l are bounded for t < a and 
(j>i{a) = ns{a, a; (f) G L;^g^(R+). In order to construct the solution of the McKendrick 
equation for any t > a, we proceed by induction and take at each step the new initial 
condition (j)i^i{a) — ns{a^a](j}i) with i > 0. This justifies the formula n{a,t) = 
ns{a, t — qa; 4>m), where q — [t/a]- D 

We introduce now a technical simplification. Defining the new function p(a, t) 
through, 

and introducing it in the McKendrick equation l|2.1(l . we obtain, 

dpja, t) ^ dpja, t) ^ ^ ^ 
dt da 

Assuming that the boundary condition for equation 14.6fl is, 

p{0,t) — / b{a,t)p{a,t) da 

J a 

with 5(a, t) = b (a, t) e" ^o Ks)ds^ ^^^^^^ ^j^^ initial condition is, 

i/,(a) = e^oX'')^«0(a) 

the Cauchy problem for the McKendrick equation (|2.1|l is simply transformed into 
the Cauchy problem for the equation 14.6|l . 

By hypothesis, for every t, 5(a, t) is zero outside the interval [a, P]. As, in general, 
the constant /3 is not an integer multiple of a, for every i, we can extend the function 
b{a,t) as a zero function in the interval [/3, /3'], where /?' = qa, g > 2 is an integer, and 
{q — l)a < /3. Hence, without loss of generality, we assume that (3 = qa, where q >2 
is an integer. In the following, and to simplify the notation, we take this approach. 

Theorem 4.2. Let n{a, 0) = 0(a) (g L;\,^(R+)) be a locally integrable initial condition 
for the McKendrick partial differential equation 1^2. Ip . with a > 0, t > 0, and boundary 
condition j^. j| ). Assume that ^(a) G Lj^^CR^) fl C^{'R-^-) is a non negative function. 
Assume further that, for every t > 0, b{a, t) is locally integrable and differentiable, 
and has compact support in the interval [a, (3], where a < (3, [3 — qa, and q > 2 is 
an integer. Define the integer m = [t/a + 1]. Then, the general solution of the Lotka 
renewal integral equation is determined recursively and is given by: 
a)IfO<t< a, 

Bi(t)= / 6(ci,i)0(ci-i)dci (4.7) 
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where Bi{t) = B{t), forO <t <a. 
b) If, a <t< 2a, 

B2{t) = b{c2,t)Bi{t-C2)dC2+ b{c2,t)^{c2-t)dC2 

J a Jt 

b{c2,t) / b{ci,t - C2)(j){ci + C2 — t)dcidc2 

J a 

+ / b{c2,t)(j){c2-t)dc2. 



c) If, ma < t < (m + l)a, with m > 2 and q > m + 1, 

Bm{t) = / b(Cm,t)Bm-l{t- Cm)dCr, 



(4.; 



rn-l pt-{m-i)a+a. 



y^ / b{Cm,t)Bm-i{t- Cm)dc„ 

.■_o J t—(m—i)a. 



f"- 

+ / b{cm,t)(t){c^-t)dc„,. (4.9) 

Moreover, B(t) is continuous, and n{a,t) is also continuous for t > (3 and a G [0,/?]. 
For t > a, the general solution of the McKendrick equation is, 

n{a,t)=TT{a)B{t-a) (4.10) 

where 7r(a), b and 4> are as in Theorem \4-.l\ 

Proof. As we have seen in the discussion preceding the theorem, it is sufficient to 
prove the Theorem for ^{a) — 0. By Theorem 14. II for m ~ 1 and t < a, we have, 

f"- 
B{t) := Bi{t) ^ b{ci,t)(j){ci-t)dci. (4.11) 



For /z 7^ 0, we make the substitution </> — > 0, as discussed above, and we obtain a). 
We consider now the case a < t < 2a. Then, 

B2{t) = / b{c2,t)n{c2,t)dc2+ b[c2,t)4i(c2 - 1) dc2 

J a Jt 

rt _ rp _ 

b{c2,t)n{0,t- C2)dc2+ / b{c2,t)(f>{c2 - t) dc2 

rt _ .p _ 

b{c2,t)Bi{t-C2)dc2+ b{c2,t)<j){c2-t)dc2 (4.12) 

where B{t) — B2{t) for a <t < 2a. Introducing 14.11|l into H4.12|l . and with (f) ^ 4>, 
we obtain b). 

With m = [t/a + 1] > 2, and assuming that /? > (m + l)a, due to the particular 
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form of the characteristic curves in Figure IXTl we obtain, 

Bm{t) = I b{Cm,t)Bjn-l{t~ c„i)dc„ 



^t — (m — 'i)Q+Q 



y^ / b{Cm,t)B,n-i(t- Cm)dCm 

j_2 Jt—(m—i)a 

n _ 

b{Cm,t)(f>{Cm - t) dCm (4-13) 



t 

where ma < t < (m + l)a and ?7i > 2, proving c). Clearly, B(t) is continuous for 
t e [0,(3], and B{t) depends on (/)(a) with a e [0,13]. As n{a,l3) = n{a)B{f3 — a), due 
to the continuity of B{t) in the interval [0, t], n{a, (3) is now continuous for a G [0, /3]. 
As B{t) remains continuous for t > (3, then n[a,t) is also continuous for t > (3 and 
ae[0,/?]. 

Defining the new initial condition ip{a) — B{(3 — a), the solution of the renewal 
equation is recursively constructed in the intervals [/3, 2(3], [2(3, 3(3], . . . 

To derive the general solution of the McKendrick equation, we use ()2.10f) . and we 
obtain the Theorem. D 

Theorems l4.1l and l4.2l give the general solutions of the McKendrick equation and of 
the Lotka renewal equation as a function of the initial data and of the boundary con- 
dition. The fertility modulus is assumed time and age dependent, and the mortality 
modulus age dependent. 

Dropping the time dependence on the fertility function b, the stability or instability 
of the solutions of the Lotka renewal integral equation and of the McKendrick equation 
follow from the previous Theorems. 

Theorem 4.3. // the fertility function is time independent and in the conditions of 
Theorem, \4-.'^ we have: 

a) If r — 1, then, for every a G [0,(3] and t > a, n{a,t) remains bounded as t —>■ oo. 

b) If r > 1, then, for every a £ [0, /3] and t > a, n{a, t) — > oo, as t —^ oo. 

c) If r < \, then, for every a £ [0,(3] and t > a, n{a,t) — > 0, as t —* oo, 

where r is the Lotka growth rate of the population, as defined in 12.11\) . In the limit 
t —^ oo, the solution of the Lotka renewal equation behaves as B{t) ~ r*''^. 

Proof. By Theorem 14.21 B{t) is continuous, and, as n{a, (3) = ■n{a)B{(3 — a), n{a,(3) 
is continuous in the closed interval [0,/3]. Therefore, there exist numbers m and M 
such that, 

m<n{a,(3)<M. (4.14) 

By Theorem 14. II for < ii < a, we have, 

n{0,ti+(3)^ B{ti+ (3)^ I b{c)n{c-ti,(3)dc. 

J a 

As, c-tiE [0, (3] for c e [a, (3] and < ii < a, by (|4.14(l . we obtain, 

mr < B{ti +(3) <Mr (4.15) 

where, by (|2.11|) . 

b{c)dc^ / b{c)e-^o^'^''>'^'dc 
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is the Lotka growth rate, and < a < /3 < cxd. Let us calculate now a bound for 
n{a,a + /?), with a € [0,/3]. By (|4.15|l . and considering the case ^ = 0, due to the 
particular form of the characteristic curves of Figure IXTl we have, 

mr < n(a, a + P) < M if r < 1 
TO < n{a,a + P) < Mr if r>l. (4.16) 



Then, by l|TT?))) . 



and, 



B{ti +a + P) ^ / b{c)n{c-ti,a + P)dc 

J a 



mr^ < B{ti+a + l3) < Mr if r<l 
mr < B{ti+a + (3) < Mr^ if r>l (4.17) 

where < ti < a. Repeating this procedure up to g = /3/a, and as, for fi = 0, 
n{a,2[3) = -6(2/3 — a), the bound for n(a, 2/3), with a G [0, /3], must obey to the 
inequalities (|4.15() . H4.17|l . etc., which gives, 

mr < n{a, 2/3) < Mr . (4.18) 

Comparing H4.14|l and (|4.18|) . the theorem follows by induction. In the asymptotic 
limit, the solution of the Lotka renewal equation behaves as B{sf3) ~ r", where s is 
an integer. D 

With 6(c) = b5{a — a) in Theorem 14.31 we obtain Theorem 13.31 
Let us see a simple example that shows that for a fertility function distributed 
along an age interval, we have damped growth cycles with a period equal to the age a. 
Suppose that 6(a) = 6 is constant in the interval [a,/? = 2a], the mortality modulus 
is age independent, and the initial population is constant in the interval (0, /3]. Let us 
denote the solution of the renewal equation by Bi{t) for t E [(i — l)a, ia], with z > 1. 
Then, by Theorem 14. 21 we have, 

Bi{t) = (/.e-'^*6a 

B2it) = (/)e~'^*6(6a(t - a) + (2a - t)) . (4.19) 

By (|^T1^ and ^U^ . we have, 

n(a,2a)=7r(a)i3(2a-a) = | ^^_,2a^^^ if a < a < 2a . ^^'^O) 

With the new initial conditions 4>{a) = n{a, 2a), by Theorem l4.2b l and b), we obtain, 
B^it) = (j)e~''%^{a^ + ^{t-2af{ba-l)) 

BJt) = c/)e~''*b'^-(bt^(ba-l)-t^(9b^a^ -6ba-3) 
6 ^ 

+ 3ta{%^a^ -ba- 8)-3a^{%^ a^ + 5ba~ 16)) . (4.21) 

In this case we have, B{0) = (t>ba, B{a) = (j)e-''°'ba, B{2a) = (l)e~'^''°'b^a'^, 
B{3a) = (j)e-^f"^b'^a^{l + 6a)/2 and B(4a) == (f)e-'^^'°'b'^a^{b + ba)/6. In Figure Ol 
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we show the time behaviour of B{t) for a = 10, b = 1, 4> = 1 and mortaUty modulus 



^ = 0.05. In this case, the Lotka growth rate is r = b{e 



-2fj,a 



)/fi = 4.77. From 



this example, we conclude that there are two time scales associated to the growth 
of new borns. The first time scale a is related with the transition from the solu- 
tions Bi{t) to Bij^iit). The second time scale (3 is associated with the Malthusian 
growth behaviour. From the above computed values of B(ia), with i > 0, we have 
B{ia) ~ r[*"/''l, in agreement with Theorem 14.31 In this example, there exists a 
damped modulation in the asymptotic time behaviour with period a. 




Figure 4.1: Time evolution of new-borns B{t) for an initial uniform population, cal- 
culated from H4.1H|I and 14.21|l . The initial condition is (/)(a) = 1, for a e (0, 2a = /3], 
and the parameter values are: 6 = 1.0, a = 10 and /i = 0.05. 

From Theorem l4.1l we can derive the Leslie discrete model of population dynamics, 
|16|. and relate the stability properties of both the continuous and discrete models. 



5 From the McKendrick model to the LesUe dis- 
crete model 

In this section, we assume that the fertility modulus is time independent. Suppose 
that the age axis is partitioned into intervals of length Aa and the extreme of these 
intervals are indexed by the integers i > 0: ai = iAa. Assume that Aa is small, in 
the sense that Aa < a. With Ai = Aa, by Theorem 14. II we have, for i > 1, 

ns {iAa, At) = ns{iAa - Aa,0)e~ Aa^- Aa A'(«)d« . 

Defining n^* = ns{iAa, At), we can write, 

n^t ^ e-,A'^^AaM(3)<i^„0_^ .^ P^-lnl, = e-AaM((»-l)Aa)„0_^ (^ > 1) (5^1) 

where we have approximated the integral by its Riemann sum, nf"^ is the density 
of individuals with age ai = iAa at time At, and Pi-i is the survival transition 
probability between age classes z — 1 and i. For i = 1, by Theorem 14. II 



.At 



ns(Aa, At) 



6 (c) (j) (c) e-'f'^" ^'^'^'^'dc . 



(5.2) 
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Approximating the integrals by Riemann sums, we obtain, 



s-l 

if* = e-^'^^^o) Aa Y^ b (mAa) (mAa) 



y^ e„(/) (mAa) 



(5.3) 



where q = [a/Aa], s = [/3/Aa], and e,; 
()5.3|l in matrix form, we obtain. 



are fertihty coefficients. Writing H5.f|) and 



,At 



,At 



/ 

Pi 








P2 



V 








es-2 Gs-l \ 







(5.4) 



Ps-2 J 



which is the discrete Leshe model for age-structured populations, [1,6' . The parametric 
relations in l|5.1|l and H5.3|l enable the comparison of both models, and in the limit 
Aa -^ 0, the solutions of the McKendrick equation converge to the solution of the 
Leslie model (|5.4|l . On the other hand, a fast way of computing numerically the 
solutions of the McKendrick equation is through H5.4|l with Aa small. 

Given an initial distribution of population numbers (n^, ■ ■ ■ , "-s-i) 7^ 0, the asymp- 
totic state of the linear system (|5.4f) is zero or goes to infinity, depending on the 
magnitude of the dominant eigenvalue of the matrix in H5.4|l . If the dominant eigen- 
value of the matrix in H5.4|l is A = 1 , bounded and non-zero population distributions 
are obtained. The eigenvalues of the matrix in (|5.4(l are determined by solving its 
characteristic polynomial equation, 0], 



P(A) = (-l)-i A-i - ^ y-'-^e, n Pj-i U . 



(5.5) 



J=2 



As we have assumed that a < f3, and Aa is sufficiently small, then e^-i > and 
es-2 > 0. Hence, the matrix in H5.4(l is primitive and, by the Frobenius-Perron 
theorem, the dominant root of the characteristic polynomial P(A) is simple, real and 
positive, |2J. Therefore, the condition for the existence of a bounded and stable 
nonzero asymptotic solution for Leslie map (|5.4|l is, 

s— 1 i 

G{Aa) := X!^* 11^^-1 ^ ^«^i ' ■■^9-1 ^ h e^.ipi . . .ps-2 = 1 (5.6) 

i=q j=2 

where the parameter G(Aa) is the inherent net reproductive number of the population 
associated with the Leslie model H5.4|l . Pj. Introducing the parameters pi and e^ as 
defined in H5.1|l and H5.3|l into 1)5. 6|) . we obtain. 



s-l 



G(Aa) = e-'^'^^WAa^6(iAa) J|, 



-AaMj-l)Aa) 



3=2 



(5.7) 



AaJ2b{iAa)e-^''^^^=" 



i;,M(0)Aa) 



i=q 
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In the limit Aa — + 0, passing from Ricmann smus to integrals, and SiS q — [a/Aa] and 
s = [/3/Aa], we obtain for the inherent net reproductive number of the population, 



limG(Aa)=/ bic)e- ^>'^''>'^'dc = r 



(5.8) 



which equals the Lotka growth rate H2.11|l . Therefore, in the limit Aa -^ 0, the Lotka 
growth rate and the inherent net reproductive number of the population are the same 
quantities. 

In the case of one fertile age class, the only non-zero fertility coefficients in H5.4II is 
Cg-i, and the matrix in (15.4(1 is non primitive but irreducible. In this case, by H5.5() . 
we have, 

P(A) = (-1)^-1 I A^-i - e._i Hp.-i J = 
and the s — 1 roots of the Leslie matrix are within the circle with radius, 

l/(s-l) 
s-1 ' '^ ' 

Therefore, when the fertility modulus is concentrated at a point, the ratio n{a, t)/NM{t) 
is an oscillatory function of time, as it has been shown in Theorem IH. 61 and NM{t) 
is the Malthusian growth function 1(3. 9() . If the fertility modulus shows dispersion 
around some age class, the matrix in 1(5.4(1 becomes primitive. In this case, for r > 1 
but close to 1, the complex eigenvalues of the matrix in 1(5.4(1 are in the interior of 
the circle with radius r and, asymptotically in time, oscillations dye out. However, 
as we have seen previously in ((4.20() , the persistence or the damping of the amplitude 
of the oscillations depend on the initial distribution of the population. 
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Figure 5.1: Fertility modulus 1(5.9(1 as a function of the dispersion parameter 7. We 
show the corresponding Lotka growth rates r calculated from 1(5. 8() , for the mortality 
modulus /x(a) = 0.001+0. 0001a, dispersion 7, and parameter values: a = 12, ao = 25, 
/? = 45 and 61 = 2. 



To study the fertility dispersive behaviour when we pass from one reproductive age 
class to several reproductive age classes, we introduce the Gaussian shaped fertility 
modulus, 



b{a) = bi 



l^-l(a- 



cof 



with 



e[a,/?] 



(5.9) 
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where 5i, aa and 7 are parameters. If 7 ^ 0, the fertiUty of the population is 
distributed among several age classes. If 7 -^ cxd, the fertility is concentrated at 
the age class a — aa, in the sense that, lini-y.^00 /_ ip{a)b(a)da — bi^(ao), where 
ip{a) € I?(R+). In Figure 101 we show the behaviour of the fertility modulus (|5.9|l 
as a function of a for several values of the dispersion parameter and consequently 
different Lotka growth rates. 

In Figure 1^?^ we show the ratios n{a, t)/NM{t) and N{t)/NM{t), where NM(t) = 
^{t~ao)/ao jg |-j^g Malthusian growth function as in (|3.9|l . but with r given by 15. 8() . 
As it is expected, for small dispersion around the class of maximal fertility, oscil- 
lations persist for a long transient time, but as we increase dispersion, their am- 
plitude decreases in time. In the limit t —^ 00, the oscillations disappear. There- 
fore, in the asymptotic limit, the solution of the McKendrick equation have the form 
n{a,t) w cir"^^*, where ci and C2 are constants and r is given by 1)5. 8|l . 
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Figure 5.2: Ratios n{a,t)/NM{t) and N{t)/NM{t) for McKendrick equation with 
boundary condition l|4.1|l . calculated from initial data (l){a) = 2, for a < 100, and 
(j){a) = 0, for a > 100. The Malthusian growth function is iVM(i) = r^t-o'o)/ao_ i^ g,) 
and c) the ratio of age distributions has been calculated at the time t — 400. a) and 
b) correspond to the fertility function of Figure [OJi), and c) and d) correspond to 
the fertility function of Figure I5.1b 'l . As it has been shown, small dispersion in the 
fertility modulus imply population oscillations. If fertility is not concentrated in one 
single age class, in the limit t — > cxd, oscillatory behaviour dies out. 
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6 Predictions from demography data 

For a population with n age classes corresponding to reproductive ages a^, with 
i = 1, . . . , n, we define the fertility function, 

n 

b{a) = ^b,S{a^a,) (6.1) 

and the boundary condition for the McKendrick equation is now, 

n 

n(0,t)=^&,n(a„i). (6.2) 

4=1 

By H2.11|l . the growth rate of the population is, 

r/3 



f p "■ "■ 

/ b{a)e- ^0 i^i^)-i^da = ^ he' ^o' ^f")'''' := ^ 
"'" i=i i=i 



r. . (6.3) 



In this case, the asymptotic behaviour of the solution of the McKendrick equation is 
determined by the Lotka growth rate of each cohort. 

Due to the linearity of the McKendrick equation H2.1|l . and the boundary condition 
(I2.2|l , the general solution is now the sum of the solutions for each r^ as in Theorem 
13.31 for time independent fertility function. In this case, each individual solution 
will have an exponential pattern of growth, modulated by a periodic function with 
period a^. Therefore, the pattern of growth of a population with several fertile age 
classes is almost periodic with several periodicities or frequencies. If the dispersion 
of the fertility around the maximal fertility age class is large, the pattern of growth 
is almost periodic in time. If, in addition, we consider time variations in the fertility 
function and mortality modulus, the pattern of growth can show strong fluctuations 
deviating from the periodic or exponential growth. This introduces a higher degree 
of unpredictability for the long time behaviour of population growth. 

One of the important aspects of these results is that the period of oscillations of 
population cycles are of the order of the age of one generation, as observed in human 
populations, |23- In the Easterlin model, the period of the cycles are of the order of 
two generations, [T^ . 

7 Conclusions 

We have obtained the weak solutions in the sense of distributions of the McKendrick 
equation and of the Lotka renewal integral equation. We have assumed an age and 
time dependent fertility modulus and an age dependent mortality modulus. 

If /x(a) is the age dependent death rate of a species, and b{a) is the age dependent 
fertility modulus, the Lotka growth rate is defined by, 

6(c)e-''o''('*)'^''dc (7.1) 

where a and /3 are the first and last reproductive age classes. The Lotka growth rate 
r determines the stability of the solutions of the McKendrick equation, in the sense 
that, asymptotically in time, we have exponential growth if r > 1, and extinction if 
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r < 1. If r = 1, we have an asymptotically stable population and the equilibrium 
solution is n{a) = npexp (— /g fJ,{s)ds), where ng is a constant. 

From the demography data point of view, (|7.1I) implies that measured fertility 
numbers by age class are given by 

b{a)e- ■f<>^'^'^' , 

partially justifying the general non symmetric shape of the fertility curves in human 
populations, ^^. 

One of the important features of the solutions of the McKendrick equation relies 
on the existence of a natural time scale determined by the age of the first fertile age 
class. In this case, a constant initial population will present cycling behaviour in the 
patterns of growth. The period of these cycles is equal to the age of the first fertile age 
class. This contrasts with the prediction of Easterlin where the period equals twice the 
mean age of one generation, ^3] . Observations in human populations corroborate this 
result, j^. These oscillations are in general damped, and the damping is proportional 
to the magnitude of the dispersion of the fertility function around the most fertile age 
class. These demography oscillations have been observed in human populations, j21| . 
and in bacterial growth in batch cultures, j2()j . 

On the other hand, it has been shown that the general solutions of the McKendrick 
equation and of the Lotka renewal equation retain the memory from initial data. 

In the limiting case of a population with only one reproductive age class, the 
modulation of the growth curve is always periodic and has no damping. The amplitude 
of oscillations increases if the mortality modulus decreases, and the limit of pure 
Malthusian growth is obtained if the mortality modulus goes to zero. The population 
retains the memory from the initial data through the amplitude of oscillations of the 
Easterlin cycles. 

It follows also from this approach that the McKendrick equation is associated 
with a small time scale, when we compare it with the time scale associated with the 
Malthusian growth law. These two time scales describe the evolution of a population 
at two different long range levels. In other words, the Malthusian growth law can 
be obtained from the McKendrick equation if the mortality is zero, the fertility is 
concentrated at one age, and the initial population is uniform along the age variable. 
Denoting by ao the age of the unique fertile age class, the population growth measured 
in time steps of ao is exponential or Malthusian. That is, the asymptotic solution of 
the McKendrick equation behaves as, 

n{a, t = kao) = n{a, ao)r('=°''-"")/"° = n{a, ao)^'^''"^^ 

for fc = 1,2,... 

We have explicitly derived the Leslie time discrete model of population dynamics 
from the general solutions of the McKendrick equation, enabling a direct calibration 
of the solutions of the McKendrick equation with population data. We have shown 
that, if the length of the age classes goes to zero, the inherent net reproductive number 
of a population, the growth parameter of the Leslie model, converges to the Lotka 
growth rate. 
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